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In a recent article, Zhang and Duan pQ presented a new potential model 
for carbon dioxide (CO2). It consists of three Lennard- Jones (LJ) sites to 
account for repulsion and dispersion and three distributed partial charges to 
describe the quadrupolar interaction. The molecular model is rigid and ro- 
tationally symmetric around the molecular axis. In that work pQ, simulation 
results on vapor-liquid equilibria (VLE), radial distribution function and self- 
diffusion coefficient have been reported for the new CO2 model which are in 
excellent agreement with experimental data. They also compared to results 
from different other models from the literature. The new model is found "to be 
superior to the previous models in general" . 

For the VLE properties, deviations between model and experimental data 
are reported to be pQ: 0.7 % for vapor pressure, 0.1 % for saturated liquid 
density, 2.3 % for saturated vapor density, and 1.9 % for heat of vaporization. 

Particularly for vapor pressure, which is the most sensitive of those proper- 
ties, the authors claim to achieve an average accuracy of better than 1 % over 
the entire temperature range with a molecular model which is noteworthy for 
any molecular model. 

'Author to whom correspondence should be adressed: Tel.: +49-711/685-66107, Fax: +49- 
711/685-66140, Email: vrabec@itt.uni-stuttgart.de 
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24 As there is widespread scientific interest in CO2, the model by Zhang and 

25 Duan pQ was employed for subsequent work [H [31 21 [5] . It should be pointed 

26 out that no VLE data were published in [51 [5J [H [51 [5] . We recently tested the 

27 CO2 model from Zhang and Duan [1] and found results that strongly deviate 

28 from these reported by Zhang and Duan lj. The new simulation results also 

29 strongly deviate from experimental data, particularly the vapor pressure and 

30 the saturated vapor density. 

31 The present assessment was made on the basis of two simulation tools that 

32 employ different methods to determine VLE. Firstly, the Grand Equilibrium 

33 method [7j was used as implemented in our simulation tool ms2 [S] and, sec- 

34 ondly, the Gibbs ensemble [S] was used as implemented in the freely available 

35 simulation tool TOWHEE [TU]. Both programs have proven to be correct e.g. 

36 in the recent Industrial Fluid Property Simulation Challenge for the case of 

37 ethylene oxide which is very similar to the present case [llj . 

38 For our approach with ms2, molecular dynamics (MD) simulations were 

39 performed in the liquid phase containing 1024 molecules. After a sufficient 

40 equilibration period, the chemical potential was calculated by Widoms's inser- 

41 tion method [T5] over 300 000 time steps. According to the Grand Equilibrium 

42 method [7j, the dew point was sampled with Monte-Carlo (MC) simulations 

43 where approximately 500 molecules were used. The simulation details were 

44 similar to those published in QT| and are not repeated here. 

45 With TOWHEE [TU], Gibbs ensemble MC simulations [§] were performed. 

46 There, smaller systems were studied, containing 500 molecules. After an equili- 

47 bration over 5 000 loops without volume and molecules transfer moves, followed 

48 by 10 000 loops with these moves, 50 000 production loops were performed. 

49 Other simulation details were similar to those published in [13] and are not 

50 repeated here. 

51 The Present simulation results are compared to those reported by Zhang 

52 and Duan 1 in Table 1 as well as in Figures 1 and 2, which also contains 

53 results from a highly accurate reference data for CO2 |14j recommended by the 
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54 National Institute of Science and Technology [15] . Figure 2 shows deviation 

55 plots where it can be seen that the present data sets based on the two different 

56 simulation methods agree with each other within their (combined) error bars 

57 throughout. However, they are significantly off the data by Zhang and Duan [I 

58 which coincide excellently with the experiment. 

59 The present simulation data (Grand Equilibrium) show the following average 

60 deviations from experimental data: 18 % for vapor pressure, 0.6 % for saturated 

61 liquid density, 17 % for saturated vapor density and 4.6 % for heat of vapor- 

62 ization. Only the saturated liquid density by Zhang and Duan [1] and the heat 

63 of vaporization are in good agreement with the experimental data, the vapor 

64 pressure and saturated vapor density are significantly too high throughout most 

65 of the temperature range. 

66 It has to be concluded that the CO2 model by Zhang and Duan [T] is not 

67 generally superior to previous models. 
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Table 1: Vapor-liquid equilibria of carbon dioxide: present simulation results with the Grand Equilibrium method (GE) [7] and the Gibbs 
ensemble (Gibbs) [9] compared to simulation results by Zhang and Duan (Zhang) pQ and the reference EOS (eos) [14]. The number in 
parentheses indicates the statistical uncertainty in the last digit. 
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Figure 1: Merker et al. 
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Figure 2: Merker et al. 
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